OUTP-05-15P 
hep-ph/0511169 



O 
O 

(N 



m 

m 
> 



o : 
Oh: 

Oh, 



X 



Nonleptonic Weak Decays of B to Dg and D mesons 

C. E. ThomafE 

Rudolf Peierls Centre for Theoretical Physics, University of Oxford, 
1 Keble Road, Oxford, 0X1 3NP 

(Dated: 23 March 2006) 

Branching ratios and polarization amplitudes for B decaying to all allowed pseudoscalar, vector, 
axial- vector, scalar and tensor combinations of Db and D mesons are calculated in the Isgur Scora 
Grinstein Wise (ISGW) quark model after assuming factorization. We find good agreement with 
other models in the literature and the limited experimental data and make predictions for as yet 
unseen decay modes. Lattice QCD results in this area are very limited. We make phenomenological 
observations on decays in to -Ds (2317) and Ds(2460) and propose tests for determining the status 
and mixings of the axial mesons. We use the same approach to calculate branching ratios and 
polarization fraction for decays in to two D type mesons. 

PACS numbers: 13.25.Hw, 12.15.Ji, 12.39.Jh 



I. INTRODUCTION 

Non-leptonic weak decays of B-mesons are impor- 
tant because they probe both electroweak physics and 
hadronic structure, and may provide a window on physics 
beyond the standard model. As well as its intrinsic in- 
terest, the hadronic part must be understood in order to 
extract electroweak and new physics from these decays. 
This involves non-perturbative QCD and so can not be 
calculated from first principles. There have been many 
attempts to model the hadronic part and their success 
varies depending on which decay mode is studied. There 
are only very limited lattice QCD results in this area. 
Semi-leptonic decays of B and D mesons have been suc- 
cessfully studied using the Isgur Scora Grinstein Wise 
(ISGW) model l]. We extend this model to the study 
of exclusive non-leptonic decays after assuming factoriza- 
tion of these decays. The model has corrections that van- 
ish when the final mesons have zero recoil and so should 
be most reliable close to this kinematic region. 

Heavy Quark Effective Theory (HQET) as used, for 
example, in Ref. Q, provides a set of symmetry rela- 
tions but not expressions for the rates. Some explicit 
model must still be used. The ISGW model satisfies the 
requirements of HQET in the zero recoil limit Q and 
provides an explicit model in which calculations can be 
made. 

There are also pole models such as that of Bauer, Stech 
& Wirbel (BSW) which we compare to our results. 
There arc different approaches for heavy to light de- 
cays such as the Soft CoUinear Effective Theory (SCET) 
framework j^] and Light-Cone Sum Rules These 
are useful where the final state mesons have high energy 
unlike the regime close to zero recoil that we are probing 
here. 

Some predictions have been made for decays to com- 
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binations of pseudoscalar Dg and D and vector D* and 
D* such as by Luo & Rosner 01 and Chen et. al. 0. 
Datta and O'Donnell j3l have made limited studies of 
the decays to axial and scalar Dg. Cheng et. al. in Ref. 
[Tll | have studied s and p-wave form factors in the light- 
front approach. They have compared this to improved 
ISGW model results given in Ref. TJ]. We specialize 
to the Ds D modes where the ISGW predictions should 
be reliable and calculate explicit results for all possible 
combinations of p and s-wave mesons consistently in one 
model giving polarization ratios where applicable. 

With B factories tightening current observations of B 
decays and finding new decay modes it is timely that 
these modes should be studied in detail. 

Our aims are to: 



• Provide robust predictions of branching ratios and 
polarization fractions for all the possible combina- 
tions. These can be compared with current and 
future experimental data. 



• Comment on what can be learnt about the nature 
of I?s(2317) and Z?s(2460) mesons and axial vector 
mixing. 



We start in SectionElwith a discussion of factorization 
and general remarks. In Section IlIII we briefly describe 
the ISGW model and extend it to non-leptonic decays af- 
ter assuming factorization. In Section Hvl we present our 
results for the branching ratios and polarization fractions 
of the allowed decay modes. We also compare our results 
with those from other models and experimental data in 
this section. Then we apply the same method to decays 
in to two D type mesons. In Section El we discuss what 
we can learn about the scalar and axial Dso and Dgi- 
We finish with some general remarks and conclusions in 
Section IVII and suggest some experimental observations 
that could be made. 
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FIG. 1: B-{hu) D-{cs)D°{cu) - Type I Tree Diagram 



II. GENERAL DEFINITIONS AND 
FACTORIZATION 

The color favored tree diagram ("Type I tree dia- 
gram") for B ^ DgD is shown in Fig. ^ Here Dg and 
D can be any s or p-wave sc and cu states. We ignore 
any contribution to this process from penguin and weak- 
annihilation topologies. 

The rate for a general Type I Tree decay B{bqi) 
Y{qiq2)X{qqi) (e.g. Y = Ds and X ~ D) can be written 
as: 

where Gp is the Fermi Constant, oi is the effective Wil- 
son coefficient, Vqb and Vq-^^q^ are Cabibbo-Kobayashi- 
Maskawa (CKM) quark-mixing matrix elements, q is the 
recoil 3-momentum in the rest frame of B and Mb is 
the B-meson mass. \A\'^ is the sum of the squares of the 
polarization amplitudes Ai: 

(2) 

i 

We use the notation i = H — , — h or II where the first 
label denotes the helicity of the X meson. 

We follow Rosner & Luo ij by assuming the factor- 
ization hypothesis (naive factorization) and writing the 
amplitude as a product of two matrix elements: 

A,^<Y\J^\Q><X\J^'\B>, (3) 

where J^ = V^ — A^ = g/7^(l — 75)9^ is the vector-axial 
current. In QCD Factorization J^] J_^J J^J it has been 
shown that B decays to two heavy mesons do not obey 
the BBNS factorization formula and so naive factoriza- 
tion is not expected to hold. However, as noted in Ref. 
|l4j . the charm quark mass is at an intermediate scale 
between the light and heavy (bottom quark mass) scales 



and so factorization is still an open question in the modes 
we consider. Empirically factorization works where it has 
been tested in these modes m 0- Comparison of our 
results with experiment can further test this assumption. 
The effective Wilson coefficient ai contains some QCD 
corrections. As remarked in Ref. this varies from 
process to process but only by less than about 1%. We 
therefore use that paper's value, ai = 1.05, for all pro- 
cesses. 

The matrix element < X\J^^\B > can be parameter- 
ized in terms of general form factors. We use the same 
definitions as ISGW and these are given in Appendix 1X1 
These form factors can then be calculated in a particular 
model such as ISGW discussed in Section ITTll The po- 
larization vectors and tensors used for J = 1 and J = 2 
mesons are given in Appendix IbI 

The matrix element < y| J^|0 > can be parameterized 
in terms of the decay constant of meson Y, /y. For vec- 
tor ('^S'l) and axial vector (^Pi, ^Pi) mesons 44] with 
polarization vector and mass My this is defined by 

<r|J^|0>=ef /yMy, (4) 

and for scalar {^Pq) and pseudoscalar (^^o) mesons 
with 4-momentum P^ , 

<y|J^|0>=zPj/y. (5) 

III. THE ISGW MODEL 

The ISGW model was developed to calculate 
semileptonic B and D decays using a constituent quark 
model and the mock meson approach p^ . We extend 
the model to non-leptonic decays using the factorization 
hypothesis and test the sensitivity to the assumptions 
made. We find that robust predictions can be made for 
the Ds D decay modes because they are in the valid kine- 
matic region. 

The matrix element < X\J^^\B > is calculated us- 
ing a non-relativistic decomposition of the quark current 
57^(1 — 75)6. This will contain corrections of 0(|q|^/m^.) 
and 0(|q|^/M^) where niq. are the quark masses appear- 
ing in the current decomposition. The model is there- 
fore reliable only close to zero recoil and with heavy 
quarks. Mb is the mass of the B meson, Mx is the 
mass of the D (or excited D) meson. My is the mass 
of the Ds (or excited Ds) meson. Mb, Mx and My 
are the respective mock meson masses. Pb, Px, Py 

are the 4-momenta, y = jrr = ^ m-^^ — ITP"; ^^^^ 

tjn = maxt = [Mb — Mx)^- The quarks are the spec- 
tator antiquark g^, the decaying quark 6, the quark in 
the final state meson with the spectator quark q, and the 
quark and antiquark pair in the other final state meson 
qi and 52- w^i, mi,, rUq, ruqi and mq2 are the constituent 
quark masses. In the decay B~ DjD^: qi — u, q — c, 
qi = s and 92 = c. 



3 



In the original ISGW paper the semileptonic differen- 
tial rate is given in terms of the form factors. We have 
adapted the model to calculate non-leptonic decays in the 
factorization hypothesis given above. The resulting rela- 
tionships between polarization amplitudes, form factors 
and decay constants are given in Appendix El For con- 
nection with the original ISGW paper [Tl, their hadronic 
tensor (Equ. 7) is given by 

V-II<^I4I^><^I^mI5>- (6) 

i 

We initially followed the original ISGW paper using 
(radial ground state) harmonic oscillator (HO) wavefunc- 
tion of the momentum space form: 

V'z=o,™,=o(k) = (7r/32)-'/'exp (-^), (7) 

V'.=i,™,=o(i)(k) = 72 (.^^)-^/^ ^exp (-^). (8) 

where (3 is the wavefunction parameter and k-^ = 
— -^(k^ + iky) . The wavefunctions are normalized such 

that / d^kij*{k)i;{k) = 1. 

The resulting ISGW form factors using harmonic os- 
cillator (HO) wavefunctions are given in Appendix ]n\ 
These are mostly the same as the ISGW paper with some 
additional ones not given in that paper. However, there 
are a couple of differences. The sign of the scalar form 
factor u+ in Equ. ID 141 is opposite to that in Equ. B37 
of Ref. As Ref. (l| does not give an expression for 
M_ it is not possible to tell if this is an overall or rela- 
tive sign. It seems likely that it is an overall sign due to 
the conventions used, in which case it has no affect on 
the results. If it is a relative sign it can only modify the 
results when meson F is a pseudoscalar or scalar. If, for 
example, Y is the pseudoscalar Dg , adding a relative sign 
changes the branching ratio significantly: from 4.2 x 10~* 
to 2.2 X 10^3 

The other difference is the additional mass ratio 
Mb/Mx ~ 2.2 in the axial vector (1^ ) form factor 
V in Equ. B43 of Ref. compared with Eau. IDlTl here. 
This is a parity violating form factor which appears in 
decays where meson y is a vector or an axial vector 
and only in the transverse polarization amplitudes. This 
has a negligible affect on the branching ratios and lon- 
gitudinal fractions. For example, when Y is the vector 
Z?*, the branching ratios (arbitrary normalization) are 
cx 5.28 X 10^'^ and 5.29 x 10^^, and the longitudinal frac- 
tions 0.930 and 0.928 respectively without and with the 
mass ratio. However, there is a significant affect on the 
parity odd fraction oc (A+_ - A^^)"^. This is 0.66 x 10~^ 
without the mass ratio but 3.1 x lO^'^ with it. 

The original ISGW paper was restricted to harmonic 
oscillator wavefunctions; our ISGW form factors for gen- 
eral wavefunctions are given in Appendix IeI 

In the original ISGW paper |qp in form factors is ap- 
proximated as {tjn — t)Mx /Mb and a relativistic correc- 



Quark Constituent Quark Mass / MeV 

b 5170 

q = qi=c 1770 

qi — u 330 
(j2 = s 550 



TABLE I: Constituent quark masses |T^ . 

tion factor k is introduced through |qp |qp/K^. This 
was intended to bring the pion form factor into better 
agreement with experiment at the low region P|. We 
take K = 0.75 18]. However, we look at how much n 
influences the results by setting k = 1 in one model. 

There is ambiguity in the mock meson approach 1^ 
as to whether, at the end of the calculation, to identify 
the mock meson masses as the physical meson masses or 
the sum of the constituent quark masses. 

We used a number of different models corresponding to 
different choices to assess the robustness of the approx- 
imations used. In model 1 the {tm — t) form is used in 
place of |qp. In models 2 and 3 the exact |qp is used 
in the exponential of F^. In models 1-2 the mock meson 
masses are equal to the physical meson masses. In model 
3 the physical meson masses are used except in the ex- 
ponential of where the sum of the constituent quark 
masses is used. Models 4-5 take the mock meson masses 
to be the sum of the constituent quark masses. Model 4 
uses {tm — t) and model 5 uses |qp. Model 6 is the same 
as model 3 except that we set n = 1. 

The Ds D decay mode is expected to be reason- 
ably well modeled by ISGW. The highest recoil momen- 
tum (that for final state pseudoscalars) is 1812 MeV. 
|q|/TO6 « 0.35, |q|/TOc « 1-02, |q|/M£, « 0.97. Although 
|q|/TOc and |q|/Af£i do not appear small, they appear 
as, for example, |q|'*/(8M|,) in a Taylor expansion with 
numerical value « 0.14. The approximation is therefore 
reasonable. In the following section we present the results 
of our calculation. 



IV. RESULTS OF CALCULATIONS 
A. ISGW Model Results 

Unless otherwise stated all numerical values of con- 
stants, masses etc. are taken from the PDG Review 2004 
|23| . The quark masses and HO wavefunction parameters 
used are shown in Tables and |H] respectively. We take 
ai = 1.05 0, \Vbc\ = 0.04 and = 0.97. The magni- 
tudes of decay constants used are given in Table UTTl The 
phases of the decay constants are chosen to match those 
given by quark model calculations in Appendix IfI 

To be explicit, we consider the mode in which a charged 
B decays in to a charged Dg and neutral . Initially 
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Meson 


13 1 MeV 


B 


410 


D 


390 


D* 


390 


Do 


330 


Di 


330 


D2 


330 



TABLE II: HO wavefunction parameters [2C|. 



Meson 


Decay Constant / / MeV 


Ds 


240 MeV [21] 


Dt 


275 MeV f21] 


Dso 


110 MeV [22j 


Dsil 


240 MeV [22] 


Dsi2 


63 MeV [22| 



TABLE III: Ds meson decay constants. 



we take the axial vector Dsi{2A60) = Dgil to be ^Pi 
with mass 2460 MeV. The 1^51(2536) = Dsi2 is ^Pi and 
has mass 2536 MeV. We discuss axial vector mixing in 
Section^ The scalar Dgo is assumed to be the Ds{2317) 
with mass 2317 MeV. We take both axial Di to have mass 
2420 MeV: Dil is ^Pi and Di2 is ^Pi. The scalar D„ has 
mass 2350 MeV and the tensor D2 has mass 2459 MeV. 

As discussed in Section IIIII to get a handle on the 
theoretical uncertainty within the ISGW model, we cal- 
culated the branching ratios and polarization ratios using 
six different model choices. 

A plot of the calculated decay rates and the experi- 
mental results is shown in Fig. |2] The model error bars 
and average show the range and average of the six model 
choices. The experimental results are from the PDG Re- 
view 2004 I2I apart from: B° DfD' which is from 



Belle "245 and B° ^ D*+D*- which is from BABAR "251. 
Tables of numerical results are given in Tables ]V\ and IVll 
of Appendix O Products of branching ratios have been 
measured for the Ds(2317) and I?s(2460) modes from 
BABAR 26] and Belle ^ . This data is averaged in the 
PDG Review 2005 Partial Update |2^ and discussed in 
Section El 

Where data exist the model seems in remarkable agree- 
ment. The possible exceptions are the newer Belle and 
BABAR results for the neutral B decaying to Dj 
and D*~ D^* . Both these neutral B modes are lower 
than the model predictions and, although this is not sig- 
nificant, they are lower than the corresponding charged 
B mode. If more precise data shows a discrepancy be- 
tween charged and neutral decays it would be due to some 
difference between u and d spectator quarks. Spectator 
interaction effects which could cause this include weak 
annihilation, which can only contribute to the charged 
5+ decay, and electromagnetic penguins. 

The variation due to the model choice is relatively 
small; it's of the same order or less than the experi- 
mental uncertainty for those modes with experimental 
data. This suggests the predictions made are robust. The 
model uncertainty shown does not include uncertainty 
in the decay constants: this would change the overall 
branching ratio for each Y (e.g. Ds) meson but not the 
pattern for a given Y . 

The Dsil D'\ Dsil D and Dsi2 D* modes are rela- 
tively large but are complicated by axial vector mixing. 
We comment on implications for Dso and Z^^i mesons in 
Section El 

Polarization ratios for the relevant decay modes are 
shown in Fig. O There is only one piece of experimental 
data in the PDG Review 2004: that from BABAR and 
CLEG on the decay in to D* D* . Here there is good 
model agreement. The model uncertainty is ~ 10% be- 
ing slightly larger for the Dil modes. Again the model 
predictions are robust. 



We then looked at the ISGW model with more real- 
istic wavefunctions using the results of Appendix IeI and 
a numerical solution of the radial Schrodinger equation. 
This was done by discretizing the position coordinate in 
to 300 intervals ranging from to w 2 x lO^^MeV^^. 
We first checked our numerical method by solving a HO 
potential V = l/2fiuj^r^ where uj = P'^/fJ. and fi is the 
reduced mass of the quarks. The results were the same 
as the HO analytic ones to the numerical accuracy. 

We then used a more realistic Linear -I- Goulomb 
-I- Hyperfine (LCHF) potential with the delta function 
smoothed out to a Gaussian 1291 [301: 



Taking = 0.594, b = 1.62 x lO^MeV^, a = 897MeV 
and an — as- The resulting form factors are neg- 
ligibly different from the HO ones compared with the 
other model uncertainties. Varying the parameters in 
the LCHF potential by 10% leads to negligible changes 
in the results. 

So, to summarize, the potential used does not appear 
to modify our results significantly. 



3 r 



32aHcr^ 



(sq • Sq) exp (-cr^r^). (9) 
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H>^ ISGW Model 

«H Experimental B° scaled by Width(B'')/Wldth(B*: 
■ <y ' Experimental 

I L_I_J I I I I I I I I I I I I I I I I I I I I I I I I I I L 



D D* Dp D,1 0^2 Dj D D* D,, D,1 D,2 Dj D D* D„ D,1 D,2 Dj D D* Dp D,1 0^2 Dj D D* D,1 0^2 Dj 

Decay Mode 

FIG. 2: B —* DsD Branching Ratios: ISGW model and experimental values. See the text for references to experimental results. 



B. Comparison with Other Results 

We compared these results with calculations based on 
the BSW model Luo and Rosner's HQET based 

model in ref. 2] and the pole model in that reference. 
We used our values of the masses, decay constants, CKM 
matrix elements and ai. These models only allow predic- 
tions when meson X is a vector or pseudoscalar meson. 
The branching fraction results are shown in Fig. ^ and 
the polarization ratios in Fig. |31 It can be seen that there 
is reasonably good agreement between the four models. 
This adds to our confidence in the robustness of our re- 
sults. The ISGW model appears to give a slightly worse 
fit in some cases. However we should note that we haven't 
fitted anything to these decays. The only parameter, ai, 
was obtained for general B decays. 

Cheng et. al. |^ have used a light-front approach 
to calculate some of these decay modes. They do not, 
however, give explicit results for all the modes that can be 
compared to experimental data. The light front approach 
has the advantage of being relativistic, however, in the 
region close to zero recoil being studied we have argued 
that the ISGW model should be valid. The light front 
approach uses a parameterization to get from maximum 
recoil to the physical regime. 



Cheng [l^l has also studied s and p-wave form fac- 
tors in an improved ISGW model Q- However different 
models are used for different transitions. For consistency 
we have used the same model for all the decay modes. 
We give explicit results for all possible combinations in- 
cluding p-wave and and p-wave modes some of which are 
comparable in magnitude to the p-wave -I- s-wave modes. 
We present polarization ratios where appropriate. None 
of the parameters have been fit to these decays; the decay 
constants are determined from other experimental results 
or models. We have also attempted to assess the impact 
of some of the approximations in our results. Their qual- 
itative pattern is the same as ours but there are quantita- 
tive differences, comparison being complicated by mixing 
effects. 

Lattice calculations currently only determine the form 
factors at zero recoil for semileptonic decays of i? to D 
or D* . Therefore branching ratios can not be compared 
with our calculations. For B ^ D the lattice results are 
given as F(w = 1) where w = Pb ■ Px/{MxMb). This 
is related to the ISGW form factor by 



^ 2^MxMb „ 

r — — , , J4 



(10) 



Mx+Mb 

The ISGW model gives 0.961 to 1.02 depending on what 
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D,2 
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D* 
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D,1 D,2 
Decay Mode 

FIG. 3: Polarization Ratios: ISGW model and experimental values. See the text for references to experimental results. 



masses are used for the mock meson masses. Lattice 
calculations give 1.074 ± 0.018 ± 0.016 [p. For B ^ D* 
the lattice results are given as Fa(w =1). This is related 
to the ISGW form factor / by 

f = ^/MxMBFA{w)il + w). (11) 

The ISGW model gives this as 0.998 whereas 
a quenched lattice calculation gives 0.913^Qg^7 ± 
0.016l:^:™^+0-noo+o.oo6 ^^^^ ^^^^^^ ^^le resuhs of 

lattice and ISGW models differ. 

C. The B ^ DD Modes 

The ISGW model can also be used in decays to two 
D (and excited) mesons. The dominant contribution is 
again expected to be from the Type I tree diagram shown 
in Fig. n In the decay B~ D~ D^: qi = u, q = c, 
qi = d and q2 — c. These decays are Cabibbo suppressed 
with the CKM matrix element Vcs replaced by Vcd leading 
to a suppression factor of 0.05 in rate compared with the 
Dg D modes. This means that other contributions may 
be relatively more important. 

An additional complication is in distinguishing the two 
D mesons produced. In neutral B decays the B^ and 



B^ decays must be distinguished to determine which D 
meson is which. In the B^ decay meson Y has charge + 
whereas in the B^ decay it has charge — . This is not a 
problem in charged B decays where meson Y is always 
charged and X is neutral. 

With these caveats the decays can be calculated as 
above. We take \Vcd\ = 0.22. The magnitudes of decay 
constants used are given in Table Hvl These are from a 
quark model calculation using HO wavefunctions, as de- 
scribed in Appendix IfI normalized to the fo result from 
CLEG •"!•> . They produce the same pattern of results as 
in Ref. 'Q. 



Meson 


Decay Constant / / MeV 


D 


223 [33] 


D* 


215 


Do 


156 


Dil 


183 


Di2 


89 



TABLE IV: D meson decay constants 

The branching ratios and polarization fractions are 
shown Figures and [T] respectively. Numerical values 
are given in Tables IVIll and IVIIII of Appendix lUl In the 
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rG -: HQET Rosner 
r Pole Rosner 

Experimental B° scaled by Wldth(B°)/Wldth(B*) 
Experimental B* 



D„1 
D 



Dj,2 
D 



D* 

Decay Mode 



FIG. 4: B ^ DsD Branching Ratios: DifTerent models. Experimental data as in Fig.|5| 
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FIG 



Decay Mode 

5: Polarization Ratios: Different models. Experimental data as in Fig. |31 



PDG Review 2005 Update [23| there are only upper Km- decay constants which could be significant. 



its for the vector and pseudoscalar modes along with a 
combined branching ratio in to _D+ D* and D of 
0.93 ± 0.15 X 10"'^. The experimental data shown on the 
figures is from Belle js^l 113 • There is good agreement 
except for decay to D D. Again here the branching 
ratio is lower than the B'^ branching ratio. The model 
uncertainty shown does not include uncertainty in the 



Because of the charge conjugation symmetry of the fi- 
nal state in decays to D* D* another interesting observ- 
able is the CP odd fraction, This is related to the 
hclicity amplitudes: Eaus. [n4l and in5l in Appendix El by 



\A+- 



2(|yl+_P + |yl_+|2 + |^„|2) 



(12) 
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H>^ ISGW Model 

«H Experimental B° scaled by Width(B'')/Wldth(B*: 
■ <y ' Experimental 

I L_!_J I I I I I I I I I I I I I I I I I I I I I I I L 



D** D*,1 D*,2 D*o 

D D* Dp D,1 0^2 Dj D D* Dg D,1 D,2 Dj D D* D,1 D,2 D D* Dp Di1 0^2 Dj D D* D,1 0^2 Dj 

Decay Mode 

FIG. 6: B ~> DD Branching Ratios: ISGW model and experimental values. See the text for references to experimental results. 



The model are shown in Table CSI of Appendix O 0.125 ± 0.044 ± 0.007 This is in general agreement 
There is experimental data on in the decay to D*D*: with = 0.086 from the ISGW model. 
Belle gives = 0.19±0.08±0.01 and BABAR i?_L = 



Chen et. al. in ref. jSj] have used generalized factor- 
ization to study decays to pseudoscalar and/or vector D 
mesons and compared a number of approaches. They in- 
clude penguin and annihilation effects and comment that 
penguin effects are not negligible in decays to two pseu- 
doscalars. They have better agreement with experiment 
in that mode and so this could be a contribution to the 
discrepancy between our branching ratio to D D and ex- 
periment. 



V. PHENOMENOLOGICAL IMPLICATIONS 

FOR Dso AND Dsi 

The L's(2317) (0+) and £13(2460) (1+) do not fit eas- 
ily in to the cs spectroscopy j37|. These states could be 
conventional cs mesons with lower than expected masses 
|38j| . Alternatively they could be multiquark or molecu- 
lar mesons associated with DK and D*K thresholds [s^ . 



Measuring branching ratios to D or D* with Ds{2317) 
or Ds{2A60) would help determine the nature of these 
mesons. If the measured branching ratios are consistent 
with our predictions this would support them being con- 
ventional mesons. If this is the case, the mixing angle of 
the "^Pi and ^Pi states in to physical axial vector mesons 
Z^s (2460) and Z^s (2536) can also be determined. 

Products of branching ratios have been measured for 
the L»j;2317) and Ds(2460) modes from BABAR '2^ and 
Belle l27.| and included in the PDG Review 2005 Partial 
Updatepjl. Following Datta and O'Donnell ^ we es- 
timate lower limits for modes containing these mesons. 

For decays to Ds(2317) D we get experimental lower 
limits of (9.0 ± 3.2) x 10"* for the B+ decay and (1.1 ± 
0.4) X 10-3 for the B° decay. For decays to 0^(2317) D* 
we obtain (9±7) x 10"* and (1.5±0.6) x lO'^ respectively. 
These are consistent with the model predictions. 

We parameterize the mixing of the axial Dg states in 
terms of a mixing angle cf). The physical states are then 
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FIG. 7: Polarization Ratios: ISGW model and experimental values. See the text for references to experimental results. 



given by: 

\Dsi) = sin0|3pi)+cos0|iPi), (13) 
iD^i) = cos^l^Fi) - sin(/.|iPi) . (14) 

When (p — 0: Dsi is purely ^Pi and corresponds to our 
1^512, and D'^^ is purely ^Pi and corresponds to our f^il. 
In the heavy quark limit the physical states are 3/^1+ 
and 1/^1+ (^3 jP j-j couphng eigenstates) corresponding 
to = -54.7°. 

Plots of branching ratios to these mesons and D and 
D* as a function of (/) are showir in Figures|Sland|51respec- 



tively. We have ignored the mass difference and assumed 
the Ds(2460) to be a conventional scalar meson. In our 
conventions the decay constants of the ^Pi and ^Pi ax- 
ials have opposite sign (see Appendix and this must 
be taken in to account when calculating the amplitudes. 
The estimated experimental limits are shown. Because 
only limits exist it is not yet possible to determine the 
mixing angle using these decays. With measured branch- 
ing ratios these decays provide a relatively clean way to 
measure the mixing angle but dependent on the D^i de- 
cay constants. 



Datta and O'Donnell (lU] claim that there is a discrep- 
ancy between experiment and theory. However they only 
consider ratios of decay constants and ignore the mass 
differences between the s and p-wave mesons. Here we 
have shown that if the ISGW model is used and the mass 
differences are taken in to account there is no inconsis- 
tency between experiment and theory. 

Ignoring effects due to the different masses between the 
■^Pi and ^Pi axial vector mesons, the only difference in 
branching ratios is due to the decay constants. The quark 



model expressions for decay constants given in Appendix 
Elimply that the branching ratio to ^Pi Dg is larger than 
that to ^Pi Dg because the first is proportional to the 
reduced mass whereas the later is proportional to the 
inverse of the difference of inverse masses. In the equal 
quark mass limit only the "^Pi meson can be produced, 
not the ^Pi. In the heavy quark limit only the ^^^1+ 
meson can be produced, not the ■^/•^1+ as remarked in 
Ref. 113 . 

The axials and vectors can not be compared to the 
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FIG. 8: Branching Ratios to Dai D as a function of mixing angle. At </!> = the top curve is D'si and the lower curve is Dsi. 
The experimental lower limits with la errors are shown for (finer dashing) and scaled _B" (coarser dashing). 




FIG. 9: Branching Ratios to Dsi D* as a function of mixing angle. At </!> = the top curve is and the lower curve is Ds\. 
The experimental lower limits with la errors are shown for (finer dashing) and scaled B^ (coarser dashing). 



scalars and pseudoscalars as simply. This is because the 
amplitudes involve different combinations of form factors 
as seen in Appendix lO 



VI. CONCLUSIONS 

The ISGW model with the factorization hypothesis 
produces robust predictions for branching ratios and po- 
larization fractions for the B — > DgD decays, which are 
in general agreement with the limited experimental data. 
The predictions are also in line with those from BSW, 
HQET and pole models. We have made predictions for 



all combinations of s and p-wave final state mesons that 
can be measured in B decays. There are opportunities 
to probe the nature of the Ds(2317) and Ds(2460) states 
and determine the axial vector meson mixing using these 
decays. Measuring more polarization fractions will be in- 
teresting because they weigh different form factors at the 
same momentum transfer. This method has also be used 
to calculate branching ratios and polarization fractions 
in B ^ DD decays. 

With the B factories at BABAR and Belle, the Teva- 
tron and later LHCb, more precise data and also data on 
other decay modes should emerge. 

It is important to note that these approximations are 
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robust for the generic class of Ds D decays but are more 
suspect for color suppressed B J/'4'K (and excited) 
decays. These are outside the kinematically valid region 
and experimental data on these modes disagree with re- 
sults from this model. In addition, the factorization hy- 
pothesis fails for these decays as evidenced by the obser- 
vation of a Xco decay mode at BABAR Ml] and Belle 
[4^ . A possible contribution to these discrepancies could 
be final state rescattering of (cs) (wc) in to (cc) {us) final 
states |43j . 
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APPENDIX A: GENERAL PARAMETERIZATION OF VECTOR-AXIAL CURRENTS BETWEEN 

MESONS 

When X is a pseudoscalar ^Sq '": 

{X \V^\ B) = f+{PB + Px)i. + f-{PB - Px)^. (Al) 

When X is a vector 1~~ with polarization vector e^: 

{X \A^\ B) = fe; + a+{e* ■ Pb){Pb + Px)^ + a.{e* ■ Pb){Pb - Px)^, (A2) 

{X \V^\ B) = ige^,p,e*%PB + PxY{Pb - PxY ■ (A3) 

When X is a scalar ^Pq 0++: 

{X \A^\ B) = u+{Pb + + u-{Pb - (A4) 

When X is an axial vector ^Pi 1++ with polarization vector e^: 

{X\V^\B)=le*^+c+{e* ■PB){PB + Px)^ + c.{e* -Pb^Pb-Px);., (A5) 

{X \A^\ B) = iqe^,p„e*%PB + PxY{Pb - PxY- (A6) 

When X is an axial vector ^Pi I"' with polarization vector e^: 

(X \V^\B) = re; + s+{e* ■ Pb){Pb + Px)^ + s.{e* ■ Pb){Pb - Px)^., (A7) 

{X \A^\ B) = ive^.p^e*''{PB + PxY{Pb - PxY, (A8) 

When X is a tensor ^P2 2++ with polarization tensor e/^u- 

{X \A^\ B) ^ ke;,P^B + b+ielpPSP^)iPB + Px)f. + 6_(e;^P^P|)(PB - Px)^, (A9) 

{X \V^\ B) = ihe^,^pe*''"PBa{PB + Px)^{Pb - PxY- (AlO) 

I 

APPENDIX B: POLARIZATION VECTORS FOR For a vector meson with 4-momentum (£, 0, 0, -|q|) 

J = 1 AND J = 2 MESONS and mass M the polarization vectors (in the spin rather 



For a vector meson with 4-momentum [E, 0, 0, |q|) and 



than the helicity basis) are 



1 



mass M the polarization vectors are = —-^(0,1,1,0), (B4) 

" (Bl) = ^{0,1,-1,0), (B5) 



" ;7f(0'l'-^'0)' (B2) = ^(-|q|,0,0,E). (B6) 



~ (|q|; 0) 0, ii^). (B3) For a tensor meson with 4-momentum (i?, 0, 0, |q|) and 
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mass M the polarization tensors can be obtained by tak- 
ing the outer product of two spin-1 polarization vectors 
with the appropriate Clebsh-Gordon coefBcients. The 
polarization tensors are given by: 



■ ■ - 





1 i 

i -1 





t-i- — 



1 

2M 





|q| 





|q| i|q| 

£; 

iE 

E iE 



(B7) 



(B8) 







r |q|^ 








|q|s 1 














^0 - 








1 

2 






1 

2 







(B9) 






L TUP" 


-0— 


-0- 


TP J 


1 



1 

2M 



|q| 

|q| 

-i|q| 



-i|q| 
E 
-iE 



E -iE 



1 

-i -1 





(BIO) 



(Bll) 



APPENDIX C: POLARIZATION AMPLITUDES IN TERMS OF FORM FACTORS 



For a pseudoscalar B decaying to a pseudoscalar (X) and a pseudoscalar (Y) the amplitude is given by 

A = ifY[{Ml-M],)f+ + My_]. 



(CI) 



For a decay in to a vector (X) and a pseudoscalar (Y) only a longitudinal vector is allowed and the amplitude is 
given by 



Al = -i^^^^ [f + a+{Ml - Ml) + a.M^] . 



(C2) 



For a decay in to a pseudoscalar (X) and a vector (Y) only a longitudinal vector is allowed and the amplitude is 

given by 



AL^-2fYf+\ci\MB. 



(C3) 



For a decay in to a vector (X) and a vector (Y) the polarizations amplitudes for X positive, negative and longitudinal 
polarization respectively are 



A+_=-fYMY{f + 2g\q\MB), 
A_+ = -/yMy if - 2fl|q|Ms) , 



(C4) 
(C5) 



All = 



Mx 



1 

4M| 



(M^ + - M^){M^ +M^- M^) + 2a+|qrM^ 



^ [fPY.Px + 2a+H''Ml] 



(C6) 
(C7) 



For a decay in to a tensor (X) and a pseudoscalar (Y) only a longitudinally polarized tensor is allowed and the 
amplitude is given by 



Al = -ifv [k + b+{Ml - M^) + b-M^] 



2A4|ql2 



3 



(C8) 



For a decay in to a tensor (X) and a vector (Y) the tensor can not have polarization ±2. The polarization amplitudes 
for X positive, negative and longitudinally polarized respectively are given by 



(C9) 
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Au 



2 HMbIy 

3 Ml 



^ (Ml + Ml ~ M^){AI^ + 



AMI 



Mi; 



2Ml\ci\H, 



(CIO) 



(Cll) 



From these expressions it can be seen that the parity violation arises from a non-zero form factors g and h. 

AmpUtudes for when Y is an axial vector or scalar can be obtained simply by using the correct decay constant. 
Amplitudes when X is an axial vector or scalar can be obtained by replacing the vector or pseudoscalar form factors 
by the analogous axial or scalar ones given in Appendix IXI anrf multiplying by —1 to account for the difference in sign 
of the vector and axial currents in V-A. 



APPENDIX D: ISGW FORM FACTORS FOR HARMONIC OSCILLATOR WAVEFUNCTIONS 



We define: 



(Dl) 



, 1 1 

/x± = — ± — 



and 



i Mx { I3bI3x 
Mb 



0" 



BX 



n/2 



exp 



tm. t 



1 4_ ™6 



The pseudoscalar (^^o) form factors are 

/+= 

f- - [l - (Mx + Mb) 

The vector {^Si) form factors are 

/ = 



1 

2mn 



4^1+ MxK 



9-^ 3 



2MbF3 

1 mg.;3| 



F3 


1 + 




2Mx 


nib 




1 - 


Mx 


2Mx 


Mb 



Mi 



The axial vector (^Pi) form factors are 



.Mb 



^5 



4t,-MBfi%x 
0%m,,Mx 
2f3l^MBl^+ 



4V2; 



= F; 



^ V2MBPE 



F. 



V2MbI3e 



I 4_ ™6 

Mx+Mb 
2mo 



nibrnqm^if}^ 



ifi+fi-MxI^l 

(Mx+MB)f3%mg, 



(D2) 

(D3) 

(D4) 
(D5) 

(D6) 
(D7) 

(D8) 

(D9) 

(DIO) 
(Dll) 
(D12) 
(D13) 



Note the difference in v compared with B43 of 0, we have checked this and believe our expression is correct. We 
comment on this in Section IlIII 



15 



The scalar (^Pq) form factors are 



rrihTaqmqi 



_ jp mgi{Mx+MB) 

5 VePsMx ■ 



Note the negative sign in u+ compared with that in equation B37 of IJ. We comment on this in Section ITTll 
The axial vector (^Pi) form factors are 



2MxPe 



^F^PbMb 



F. 



-F 



m,i{MB+Mx) 
4Mb 13 b 



1 rrigiPj 



The tensor (^P2) form factors are 

h = 
k = 

b- 



F. 



2^MbI3e 



V2F,^, 



^ 2V2Mxmbl3E 

F. 



2,i+MbP%x 



1 



2y/2MlPB 



-1 



Mb 
Mx 



iMBli-l3%x V -^MbP^bx 
_ niqi , mgi{Mx+MB) 



(D14) 
(D15) 



(D16) 
(D17) 
(D18) 
(D19) 

(D20) 
(D21) 
(D22) 
(D23) 



APPENDIX E: ISGW FORM FACTORS FOR GENERAL WAVEFUNCTIONS 



We define: 



and 



/J|xs j(/5|+«). 



nib 



(El) 



(E2) 



In the following () = {X\B) and (O) = {X \0\ B) . For example (p^) = J d'^p.pz.tp* {p)ip{p). For notational convenience 



we include the common y Mx / Mb factors in () and (O). 
The pseudoscalar {^Sq) form factors are 



f- 



(i + i^) 
Ofi- 



Mx+Me 
2m o 



jpz) rn,,mg 
|q| 2/i_/j+' 

(Pz) Mx+Mb 



|q| 



2m+ 



The vector (^5*1) form factors are 

9 = 
f = 







nig |q|/i_ 

2(>Mb, 



/\ ^^ -i-Mjl 



Mi 



(E3) 
(E4) 



(E5) 
(E6) 
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The axial vector ("'^Pi) form factors are 







|q|A/B 



nibmqMx _ (]5+) 



|q|M£ 



1 - 



2m o 



2A^BM+|qP 



1^ 



2^_Afx 

(A/x+Mb)A?x 
2A/BA'+|q|2 • 



(E9) 
(ElO) 

(Ell) 
(E12) 



The scalar (^Pq) fo rm factors are 



_ A 1 rating 

- V 73 |q|M~ ' 

_ /\ 1 (Mx+Mb) 

^ ^' Vs. Iql 



(E13) 
(E14) 



The axial vector ('^Pi) form factors are 



c+ = 



N/2^|q| ' 



Mb 



(p.) I ()|q| I (P+) 
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2V2A/B|qP 
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(E15) 
(E16) 
(E17) 
(E18) 



The tensor (^^2) form factors are 

k = 
h = 

h+ = 



Q 2Afx 



A-fx 



2A/B|qP 



|q| ' 
, (>|q| (p+ 



Mi 



2|qPAf| 



Mi 



2|q|2M^ 



2/^_ 



2/i_ 
(Mx+Mb) 



mi,mq 
Af2 ^_ 
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(E19) 
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APPENDIX F: DECAY CONSTANTS IN THE 
QUARK MODEL 

A non-relativistic quark model calculation using simple 
harmonic oscillator wavefunctions gives: 



hp. 
hp, 

hp. 
hp. 



= i\/2' 



N{Y) 



V2. 



My 
N{Y) 
My ' 



PyN{Y) 



H^My 
_^PyN(Y) 
^1+My ' 
I3yN{Y) 
^i^My ' 
0. 



(Fl) 

(F2) 

(F3) 

(F4) 

(F5) 
(F6) 



/3y is the SHO wavefunction parameter, 
N{Y) cx (47r/3^)3/My 

and 

± = ^±^. 

fl± mq2 rUql 

Here to^i is the outgoing quark mass and 771^2 is the 
outgoing antiquark mass. The sign of the ^Pi decay con- 
stant depends on which way around the spin and orbital 
angular momentum are coupled. We use the convention 
consistent with the argument in Ref. this work, a 

different relative sign here would only lead to the branch- 
ing ratios shown in Figs. IHlfcElbeing shifted in 

For general wavefunctions the above equations are re- 
placed by: 
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f^Po 
/a Pi 
f^P. 

hp. 



V2^^^?^ < 0,0|0 >L, 



My 

1 Ar(y) 



1 iv(y) 



(< i,o|p,|o >L + < 1, i|p+|o >L + < 1, -ib-|o >l) , 



2 ^+Afy 

1 iv(r) 

\/2/i-My 

0. 



(<l,lb+|0>L + <l,-lb-|0>L): 
< l,0|p,|0 >L, 



Here the overlaps < l,mi\pi\Q >l only contain the spa- 
tial integrals; the spin and Clebsh-Gordan factors have 
already been taken care of. p+ = —^[p,j; + ipy) and 

P- = -^{Px-iPy)- 

Note that a tensor meson can not be produced from 
an axial- vector current. 

In the equal mass limit ^ = and so the scalar or 
^Pi axial can not be produced. However, the ^Pi axial 
can be produced. 

In the heavy quark limit where the quark mass m^i — > 
oo and p^ = fj.+ both axial vectors can be produced. It 
is useful to change from the L-S basis '^"'^^Lj to the j-j 
coupling basis ^ where j is the total angular momen- 
tum of the light degrees of freedom, J is the total angular 
momentum and P is the parity. The transformation is 
given by j4Qi |: 



\'Pi) 



1/2-^ 



3/2 + 
3/2j^ 



(F13) 
(F14) 



(F7) 

(F8) 

(F9) 

(FIO) 

(Fll) 
(F12) 



r 



1/2^- 



^rA)-V^i-M) (F16) 



In the heavy quark limit /.^p^ — — V2/ipi and so: 



/1/20+ 



0, 



(F17) 
(F18) 
(F19) 



So the 3/2j^+ axial can not be produced but the 
axial can. Note also that, apart from a phase, the scalar 
has the same decay constant. This is a specific realization 
of the general result given by Le Yaouanc et. al. in Ref. 

m- 

It is important to note that if it was the antiquark mass 
that was large there would be a relative negative sign be- 
tween /i_ and /i+ which would change the negative signs 
around in the basis change between L-S and j-j coupling. 
This means that the basis transformation has different 
signs in Df and Dj mesons. 



J- J- APPENDIX G: TABLES OF RESULTS 



Decay Mode 


Branching Ratio 


Lower 


Upper 


d: d 


1.49 X 10" 


2 


1.27 X 10" 


2 


1.70 X 10" 


-2 




3.18 X 10" 


2 


2.84 X 10" 


2 


3.50 X 10" 


-2 


D: Do 


6.97 X 10" 


4 


6.41 X 10" 


4 


7.36 X 10" 


-4 


D* Dil 


1.54 X 10" 


3 


1.20 X 10" 


3 


1.70 X 10" 


-3 


D* Di2 


2.20 X 10" 


3 


1.79 X 10" 


3 


2.51 X 10" 


-3 


Dt D2 


1.34 X 10" 


3 


1.06 X 10" 


3 


1.54 X 10" 


-3 


Dsil D 


8.10 X 10" 


3 


7.18 X 10" 


3 


8.92 X 10" 


-3 


Dsil D* 


2.67 X 10" 


2 


2.45 X 10" 


2 


2.88 X 10" 


-2 


Dsil Do 


2.70 X 10" 


4 


2.55 X 10" 


4 


2.84 X 10" 


-4 


D,il Dil 


6.82 X 10" 


4 


5.68 X 10" 


4 


7.25 X 10" 


-4 


Dsil Di2 


7.17 X 10" 


4 


6.13 X 10" 


4 


7.88 X 10" 


-4 


Dsil D2 


5.07 X 10" 


4 


4.21 X 10" 


4 


5.61 X 10" 


-4 


Dsi2 D 


5.07 X 10" 


4 


4.53 X 10" 


4 


5.53 X 10" 


-4 


Dsi2 D* 


1.86 X 10" 


3 


1.72 X 10" 


3 


2.00 X 10" 


-3 


Dsi2 Do 


1.48 X 10" 


5 


1.40 X 10" 


5 


1.55 X 10" 


-5 


Dsi2 Dil 


3.79 X 10" 


5 


3.21 X 10" 


5 


4.00 X 10" 


-5 


Dsi2 Di2 


3.84 X 10" 


5 


3.33 X 10" 


■' 


4.18 X 10" 


-5 


Dsi2 D2 


2.61 X 10" 


5 


2.19 X 10" 


5 


2.87 X 10" 


-5 


Ds D 


1.76 X 10" 


2 


1.49 X 10" 


2 


2.05 X 10" 


-2 


Ds D* 


7.51 X 10" 


3 


6.38 X 10" 


3 


8.88 X 10" 


-3 


D, Do 


4.63 X 10" 


4 


4.21 X 10" 


4 


4.94 X 10" 


-4 


Ds Dil 


2.16 X 10" 


4 


7.77 X 10" 


5 


2.67 X 10" 


-4 


Ds Di2 


2.76 X 10" 


3 


2.19 X 10" 


3 


3.20 X 10" 


-3 


Ds D2 


3.21 X 10" 


4 


2.15 X 10" 


4 


3.92 X 10" 


-4 


Dso D 


3.36 X 10" 


3 


2.97 X 10" 


3 


3.81 X 10" 


-3 


Dso D* 


1.28 X 10" 


3 


1.12 X 10" 


3 


1.50 X 10" 


-3 


Dso Do 


4.98 X 10" 


5 


4.64 X 10" 


5 


5.27 X 10" 


-5 


Dso Dil 


1.57 X 10" 


5 


7.04 X 10" 


6 


1.87 X 10" 


-5 


Dso Di2 


3.05 X 10" 


4 


2.49 X 10" 


4 


3.44 X 10" 


-4 


Dso D2 


2.44 X 10" 


5 


1.69 X 10" 




2.91 X 10" 


-5 



TABLE Y: B ^ DsD branching ratios calculated in the 
ISGW model. The branching ratio, lower and upper values 
give the average and range of results dependent on model 
choices as discussed in Section UTTl 



Decay Mode 


Longitudinal Polarization Fraction 


Lower 


Upper 


Dl D* 


0.47 


0.46 


0.49 


Dt Dil 


0.19 


0.10 


0.22 


D* Di2 


0.93 


0.92 


0.93 


D: D2 


0.50 


0.48 


0.50 




0.42 


0.42 


0.44 


Dsil Dil 


0.15 


0.10 


0.16 


D,il Di2 


0.81 


0.80 


0.82 


Dsil D2 


0.44 


0.44 


0.45 


Dsi2 D* 


0.41 


0.41 


0.43 


D,i2 Dil 


0.15 


0.11 


0.16 


Dsi2 _Di2 


0.76 


0.75 


0.77 


Dsi2 D2 


0.43 


0.43 


0.44 



TABLE VI: B -* DsD polarization fractions calculated in 
the ISGW model. The polarization fraction, lower and upper 
values give the average and range of results dependent on 
model choices as discussed in Section ITTTl 



Decay Mode 


Branching Ratio 


Lower 


Upper 


D*+ D 


5.04 X 10" 


4 


4.26 X 10" 


4 


5.81 X 10" 


-4 


D'+ D* 


9.58 X 10" 


4 


8.48 X 10" 


4 


1.06 X 10" 


-3 


D*+ Do 


2.51 X 10" 


5 


2.29 X 10" 


5 


2.67 X 10" 


'5 


D*+ _Dil 


5.29 X 10" 


5 


4.00 X 10" 


5 


5.90 X 10" 


'5 


D*+ Di2 


8.35 X 10" 


5 


6.71 X 10" 


5 


9.61 X 10" 


'5 


D*+ Da 


4.70 X 10" 


5 


3.66 X 10" 


5 


5.45 X 10" 


-5 


D+l D 


2.54 X 10" 


4 


2.24 X 10" 


4 


2.81 X 10" 


-4 


D+l D* 


7.93 X 10" 


4 


7.27 X 10" 


4 


8.58 X 10" 


-4 


D+l Do 


8.98 X 10" 


6 


8.45 X 10" 


6 


9.43 X 10" 


-6 


D+l Dil 


2.24 X 10" 


5 


1.85 X 10" 


5 


2.39 X 10" 


-5 


D+l Di2 


2.42 X 10" 


5 


2.06 X 10" 


5 


2.67 X 10" 


-5 


D+l Da 


1.71 X 10" 


5 


1.42 X 10" 


5 


1.91 X 10" 


-5 


D+2 D 


6.00 X 10" 


5 


5.30 X 10" 


5 


6.64 X 10" 


'5 


D+2 D* 


1.88 X 10" 


4 


1.72 X 10" 


4 


2.03 X 10" 


-4 


D+2 Do 


2.12 X 10" 


6 


2.00 X 10" 


6 


2.23 X 10" 


-6 


D+2 Dil 


5.30 X 10" 


6 


4.38 X 10" 


6 


5.66 X 10" 


-6 


D+2 Di2 


5.72 X 10" 


6 


4.86 X 10" 


6 


6.32 X 10" 


'6 


D+2 Da 


4.05 X 10" 


6 


3.35 X 10" 


6 


4.51 X 10" 


-6 


D+ D 


7.95 X 10" 


4 


6.69 X 10" 


4 


9.36 X 10" 


-4 


D+ D* 


3.46 X 10" 


4 


2.92 X 10" 


4 


4.12 X 10" 


-4 


D+ Do 


2.36 X 10" 


5 


2.13 X 10" 


5 


2.54 X 10" 


-5 


D+ Dil 


1.19 X 10" 


5 


4.10 X 10" 


6 


1.48 X 10" 


'5 


D+ Di2 


1.39 X 10" 


4 


1.09 X 10" 


4 


1.62 X 10" 


-4 


D+ Da 


1.72 X 10" 


5 


1.14 X 10" 


5 


2.12 X 10" 


-5 


D+ D 


3.44 X 10" 


4 


3.05 X 10" 


4 


3.88 X 10" 


-4 


D+ D* 


1.29 X 10" 


4 


1.13 X 10" 


4 


1.50 X 10" 


-4 


D+ Do 


4.73 X 10" 


6 


4.42 X 10" 


6 


5.01 X 10" 


'6 


D+ Dil 


1.42 X 10" 


6 


6.55 X 10" 


7 


1.68 X 10" 


-6 


Do+ Di2 


2.89 X 10" 


5 


2.37 X 10" 


5 


3.25 X 10" 


-5 


Do+ Da 


2.19 X 10" 


6 


1.52 X 10" 


6 


2.60 X 10" 


-6 



TABLE VII: B DD branching ratios calculated in the 
ISGW model. The branching ratio, lower and upper values 
give the average and range of results dependent on model 
choices as discussed in Section UTTl 



Decay Mode 


Longitudinal Polarization Fraction 


Lower 


Upper 


D*+ D* 


0.49 


0.48 


0.51 


D*+ Dil 


0.22 


0.11 


0.24 


D*+ Di2 


0.94 


0.94 


0.95 


D*+ D2 


0.51 


0.50 


0.52 


D+1 D* 


0.43 


0.42 


0.44 


D+1 Dil 


0.15 


0.098 


0.17 


D+1 Di2 


0.83 


0.82 


0.84 


D+1 D2 


0.45 


0.44 


0.45 


D+2 D* 


0.43 


0.42 


0.44 


D+2 Dil 


0.15 


0.098 


0.17 


D+2 Di2 


0.83 


0.82 


0.84 


D+2 D2 


0.45 


0.44 


0.45 



TABLE VIIL B DD polarization fractions calculated in 
the ISGW model. The polarization fraction, lower and upper 
values give the average and range of results dependent on 
model choices as discussed in Section ITTTl 



Decay Mode 


CP Odd Fraction, R±_ 


Lower 


Upper 


D*+ D* 
D*+ Dil 
D*+ Di2 
D'+ Da 


8.61 X 10"^ 
4.55 X 10"^ 
5.76 X 10-* 
5.82 X 10"^ 


7.93 X 10"^ 
3.99 X 10"^ 
5.69 X 10"* 
5.68 X 10"^ 


8.95 X 10"^ 
6.15 X 10"^ 

5.89 X 10"* 

5.90 X 10"^ 


D^l D* 
D+1 Dil 
D+1 Di2 
D+1 Da 


7.12 X 10"^ 
5.00 X 10"^ 
9.33 X 10"* 
3.62 X 10"^ 


6.65 X 10"^ 
4.56 X 10"^ 
9.32 X 10"* 
3.45 X 10"^ 


7.36 X 10"^ 
6.05 X 10"^ 
9.35 X 10"* 
3.71 X 10"^ 


D+2 D* 
D+2 Dil 
D+2 Di2 
D+2 Da 


7.12 X 10"^ 
5.00 X 10"^ 
9.33 X IQ-" 
3.62 X 10"^ 


6.65 X 10"^ 
4.56 X 10"^ 
9.32 X 10"* 
3.45 X 10"^ 


7.36 X 10"^ 
6.05 X 10"^ 
9.35 X 10"* 
3.71 X 10"^ 



TABLE IX: B DD CP odd fractions Ri_ calculated in 
the ISGW model. The polarization fraction, lower and upper 
values give the average and range of results dependent on 
model choices as discussed in Section ITTTl 



